Quantum non-triviality in two dimensional p-orbital Bose-Einstein condensates 
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We investigate the quantum phase transition in a p-orbital bosonic system using Wilsonian renor- 
malization group, where the p-orbital bosons condense at nonzero momenta and display rich phases 
including both time-reversal invariant and broken BEC states. We show that a quantum phase tran- 
sition between the two different condensed states can be induced by the quantum fluctuations of 
the interactions between p x and p y orbital bosons even if the interactions are marginally irrelevant. 
This is a phenomenon beyond the well-known result of "quantum triviality" which states that in two 
dimensional boson systems the boson-boson repulsive interactions are marginally irrelevant so that 
the quantum critical behaviors of the systems are noninteracting, thus "trivial". This non-trivial 
phenomenon can be explained as an effect of the quantum fluctuation induced symmetry breaking. 
We also discuss the observation of these phenomena in the realistic experiment. 

PACS numbers: 03.75. Nt, 05.30. Jp, 67.85. Jk, 67.85.Hj 



Introduction. — Confining cold atoms in an optical lat- 
tice has proven to be an exciting and rich environment 
for studying many areas of physics [lH5[- However, the 
ground state wavefunctions of single component bosons 
are positive defined in the absence of rotation as de- 
scribed in the "no-node" theorem [6J, which imposes 
strong constraint for the feasibility of using boson ground 
states to simulate many-body physics of interest. One 
way of circumventing this restriction is to consider the 
high orbital bands since the "no-node" theorem only ap- 
plies to ground states The unconventional Bose- 
Einstein condensations (BECs) of high orbital bosons ex- 
hibit more intriguing properties than the ordinary BECs, 
including the nematic superfluidity [1, 0], orbital su- 
perfluidity with spontaneous time-reversal (TR) symme- 
try breaking 1(3- 1(| and other exotic properties frH - fljj ]. 
The theoretic work on the p-orbital fermions is also ex- 
citing |20h25| . Furthermore, the p-orbital and multi- 
orbital superfluidity have been recently realized exper- 
imentally by pumping atoms into high orbital bands [j§- 
SEElIa- Hemmerich's group has successfully populated 
87 Rb atoms into p-orbital bands in a double-well optical 
lattice with a life time of several tens of milliseconds 0] ■ 
The condensate momenta are located at half values of 
the reciprocal lattice vectors, which is consistent with 
the appearance of staggered ordering of orbital angular 
momenta. 

Since the p-orbital Bose gas exhibits rich phase struc- 
tures it is interesting to explore the quantum phase tran- 
sition in such a system. In a general two dimensional 
boson system the problem of Quantum triviality exists. 
Quantum triviality was first discussed by L. D. Landau 
28] . In quantum electrodynamics the perturbative renor- 
malization group (RG) calculation shows that the renor- 
malized charge vanishes as the cutoff sent to infinity. The 
theory is noninteracting, and thus "trivial" . This indi- 
cates the failure of the perturbative theory. K. Wilson 
and J. Kogut also suggested that the c/> 4 quantum field 



theory is in fact trivial in dimensions d > 4 [29( . Espe- 
cially, at the upper-critical dimension d — 4 the triviality 
depends on whether the interaction is repulsive or at- 
tractive. If it's repulsive interaction, the renormalized 
coupling goes to zero as the length scale approaches in- 
finity, therefore the theory is trivial. A paradigmatic de- 
scription of dilute Bose gas is a <p A theory with dynamic 
exponent z — 2. At zero temperature the upper-critical 
dimension of such a system is d = 2. The renormalization 
group analysis shows that the boson-boson interaction is 
marginally irrelevant at one-loop level in two dimensions 
if the interaction is repulsive [3fJ-|32| . Hence, the s-orbital 
boson system is bound to quantum triviality in two di- 
mensions. The p-orbital boson system is beyond the "no- 
node" paradigm, whose ground state wavefunctions can 
be complex- valued. It has richer interaction properties 
and quantum fluctuation effects than the s-orbital one, 
and hence may exhibit nontrivial behaviors beyond the 
quantum triviality. 

In this Letter, we demonstrate that quantum trivial- 
ity is genuinely nontrivial for a two dimensional p-orbital 
bosonic system in a square lattice. It turns out the inter- 
actions between the p x and p y orbital bosons can generate 
nontrivial effects even if they are marginally irrelevant. 
Before they renormalize to zero, these interactions can 
modify the flow of the chemical potentials and eventu- 
ally change their flow directions, which lead to a quan- 
tum phase transition of the system. The spectrum of p- 
orbital bosons in square lattice has two energy minima in 
the Brillouin zone located at K\ — (~, 0) for p^-band and 
K% = (0, -) for p a -band respectively @, 0- At these two 
band minima the Bloch wave functions are time-reversal 
invariant and, thus, real valued. Lattice asymmetry fa- 
vors a ground state at either K\ or K%, which is called 
real BECs. A linear superposition of these two real val- 
ued wave functions with a fixed phase difference forms a 
complex BEC, which is favored by the system with inter- 
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species interactions between p x and p y orbital bosons and 
spontaneously breaks the time- reversal symmetry 0, 1X3 . 
In our work we follow the spirit of "Shankar's RG" [33] 
to investigate if any quantum phase transition between 
the real and complex BEC phases can be induced when 
various interactions are turned on. Before proceeding to 
detailed calculation, we give a simple physical explana- 
tion of our result, which is illustrated in Fig. (p}. We 
find that without inter-species interactions the complex 
BEC phase is confined in the first quadrant of the phase 
diagram as shown in (a) of Fig. (fTJ). This represents 
a non-interacting system. Quantum triviality survives. 
However, as we turn on the inter-species interactions the 
complex phase is enlarged into the second and fourth 
quadrants as shown in (b), (c) and (d) of Fig. (p}. That 
is, certain region of real BEC phases become unstable 
when the inter-species interactions are turned on and the 
system finally flows to the complex BEC phase. These 
interactions give rise to an instability from the real BEC 
phase to the complex BEC phase. 
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FIG. 1. (Color online) Phase diagrams of the p-orbital boson 
system. Four phases are determined by the flow directions 
of the chemical potentials [i\ and p,2 of p x and p y orbital 
bosons under the renormalization group transformation. In 
the complex BEC phase p,i(£) and fti{t) run to positive in- 
finities, which represents that both p x and p y orbital bosons 
condense. In the real BEC 1 (or 2) phase the p x (or p y ) or- 
bital boson condenses and the other one doesn't. In the no 
BEC phase neither of the two orbital bosons condense. The 
inter-species interaction coupling 33(0) between the p x and p y 
orbital bosons is set to 0,0.1,0.3 and 0.5 for graphs (a), (b), 
(c) and (d) respectively, "bl" and "b2" denote boundary 1 
and boundary 2 between the real and complex BEC phases. 

The renormalization group flow equations and phase 
diagrams. — The tight-binding model of the p-orbital 



bosons in a square lattice shows that the energy minima 
are located at half values of the reciprocal lattice vectors 
d, 0). Around these two band minima K\ = (— ,0) and 
K2 = (0, — ) the system can be approximated as an inter- 
acting two-component Bose gas. The action of this low 
energy effective theory with the most general interactions 
is cast as following 
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(1) 



= JZo t Jo ^« fe )H" + ^ *) is 

the dynamic term of the boson field 4>i(i = 1,2). Here 
we give a cutoff A to the momentum space since this is 
a low energy effective theory, fii and [12 are chemical 
potentials of the low energy modes <j>\ and <p2- The 
self-interaction term of each boson field is written as 
Si = J ijk g l (f>*(uj4,k4 : )(j)*(uj 3 ,k 3 )(p t (uj2,k2)(p t {uJi,k 1 ). 

Here we used a short-handed notation J , — 

Ilti Ho ^ Jo &?(2K) 2 S(k 4 +k 3 -b 2 ~kM^)S(u 4 + 
uj 3 — L02 — oji). The most general inter-species interac- 
tions between <j>i and 4>2 include two types of couplings 
S i2 = fuk {.93^2(^4, ki)4>*{u$, k 3 )4>i(u>2, k 2 )(f>2(^i,ki) + 

S>4[0i(w4,fei)0i(w3, £3) 02(k>2, k 2 ) 02 (wi,fci) + H.C.]\. 
The c/4 term can rise in this high orbital model as 
an Umklapp scattering process since the momentum 
transfer is ±2(Ki — K 2 ), which equals to the reciprocal 
lattice vectors 0, (3) ■ 

In contrast to the mean-field theory where certain 
order parameter is presumably defined, the renormal- 
ization group analysis treats various instabilities on an 
equal footing without assuming any specific order param- 
eters. We start from a microscopic theory and study the 
running of the couplings and chemical potentials using 
momentum-shell renormalization group method. Follow- 
ing the Wilson's approach [29( one finds that the chemical 
potentials scale according to the following relations up to 
one-loop order: 

^ = 2£i + 4& • 0(£i - 1/2) + ~g 3 ■ 0(£ 2 - 1/2), 

^ = 2^2 + 45 2 • e{fi 2 - 1/2) + 33 • 9{fix - 1/2). (2) 

In above equations we defined the dimensionless chem- 
ical potentials and coupling constants as jli — fam/A 
and gi = <7,m/(27r). The 6 functions in the chemi- 
cal potential flow equations are from uj integrations in 
the one loop calculations, where we have 6(p, — 1/2) = 



+00 dui e*" 
■00 2tt iuj — (1/2— /i) 



Notice that we introduced the factor 



into the ui integral otherwise the integral over u 



doesn't converge [33j. Most of the studies investigated 
the behaviors around the critical point, then these con- 
tributions can be ignored at zero temperature 31 , [32| . 
However, in our case we focus on the asymptotic behav- 
iors of the chemical potentials. The running behaviors of 
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Hi and fi2 to +00 or —00 are used as the criteria of which 
phase the system will fall into. Therefore, we have to take 
these contributions into account. They will give severe 
influence to the running of the chemical potentials. The 
flow equations of the coupling constants are as following: 



d£ ~ ^ ' 

dg3 

de 



-ah 



dg 2 

de 
de 



-252 - 2gi 

-2<?i<?4 - 2g 2 g 4 . (3) 



All the coupling constants with positive initial values are 
marginally irrelevant. For example, 173 can be solved as 



93(0) 



It approaches to zero as the length 



1+93 (o)r 

scale e goes to infinity. However, it doesn't imply that 
we can ignore these irrelevant coupling constants. This 
is because the small gi will generate small contributions 
to jEtj, which will then quickly grow under the renormal- 
ization. As discussed by R. Shankar 33j, an irrelevant 
operator can modify the flow of the relevant couplings 
before it renormalizes to zero. 

The running chemical potentials fii(£) and jl 2 {£) are 
relevant and can be solved numerically by Eq.([2]). Wc 
find that in regions of /2i(0) > \ & ^2(0) > \ and 
Ai(0) < & /i 2 (0) < 0, the running directions of p,\{t) 
and jl 2 {£) are completely determined by their initial val- 
ues. If the initial values are positive or negative, they 
finally flow to positive or negative infinity. However, in 
other regions the running directions of the chemical po- 
tentials can eventually be changed by the one-loop correc- 
tions from the interaction couplings, even if they renor- 
malize to zero. For instance, in Fig. ([2]) we start the 
running of fM\(£) from a negative initial value. As we 
vary the interaction coupling 33(0) from 0.3 to 0.8 we 
observe that the running of /2i (£) can finally be changed 
from the negative to the positive direction. That is, even 
if fix (£) runs to negative infinity at the tree level, the pos- 
itive one-loop contributions can make jli (£) go to positive 
infinity eventually. The system will finally end up in a 
different phase. 

Based on the numerical calculations the phase dia- 
grams can be drawn in Fig. ([1]). The four phases are 
determined by the flow directions of fli(£) and jl 2 (e) as 
e — > 00. (I) Complex BEC: jli(£) — > +00 and fail) — > 
+00, (II) Real BEC 1 : fn{£) -)• +00 and fa{t) -» -00, 

(III) Real BEC 2: £n(£) -)■ -00 and £l 2 (£) -> +00 and 

(IV) No BEC: jli(£) -> -00 and fi 2 (£) -> -00. In Fig. 
(fT]) the initial values of the couplings are <?i (0) = g 2 (0) = 
34(0) = 0.1 for all of the four diagrams. 33(0) is given in 
the caption of Fig. (fT]) . Comparing the four phase dia- 
grams in Fig. (fT]), we see that as the inter-species interac- 
tion 173(0) becomes stronger the complex BEC phase get 
enhanced. Given that 173(0) is small, we can obtain the 
approximate expressions of the two boundaries. Bound- 



ary 1 is /I 2 (0) 
and j&2(0) = 
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FIG. 2. (Color online) The flow directions of chemical poten- 
tial p,i(£) with different interaction coupling 33(0). As 33(C)) 
is increased the running direction of jii (£) changes from nega- 
tive infinity to positive infinity. The system can finally flow to 
a condensed phase. The initial values of the parameters are 
£i(0) = -0.125, /i 2 (0) = 0.3 and <?i(0) = g 2 {0) = g 4 {0) = 0.1. 
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> 



for < j5 2 (0) < ^ and 
j. They are indicated by 



Mi(0) = -MS for jl 2 (0) 
"61" and "62" in Fig. ©. 

Quantum fluctuation induced symmetry breaking 
( QFISB). — The mean-field results of the phase transition 
was derived by constructing a Ginzburg-Landau theory 
in Ref. [141 ] . However, starting from a microscopic theory 
our renormalization group analysis gives some qualitative 
differences: (I) In the mean-field analysis the boundary 
conditions between the real and complex BEC phases 
depend on the self- interaction couplings g\ and g 2 ■ How- 
ever, these two couplings don't affect the phase bound- 
aries in our results. The quantum phase transition is 
purely induced by the inter-species interaction between 
the p x and p y orbital bosons in RG analysis. (II) At 
one-loop level the 174 term doesn't give any contributions 
to the flow equations of the chemical potential \i\ and 
/i 2 . It can get involved in higher order calculations. For 
instance, 34 term can generate corrections to the bound- 
aries "bl" and "b2" at two-loop level through the sunrise 
graphs in Fig. ([3]). However, in mean-field analysis the 
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FIG. 3. The Feynman diagrams of the two-loop corrections 
to the chemical potential (a) Hi and (b) fi2- <j>i and 
the boson fields defined in Eq. (1). 



t>2 are 



phase boundaries depend on both 174 and 33. This differ- 
ence originates from the starting points of the two analy- 
sis. The 34 term explicitly breaks the original U (1) x U(l) 
symmetry to Ud(1) symmetry, where the index "D" de- 
notes for "Diagonal" , and leads to a fixed phase difference 
between the two fields <pi and 4> 2 . The mean-field analy- 
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sis starts from this symmetry breaking phase. Hence, the 
complex phase in mean-field is a coherent superposition 
of the two ground states. However, our renormalization 
group analysis starts from the normal phase and focuses 
on the effects of the quantum fluctuations. In this case g± 
term doesn't show its contributions up to one-loop level. 
Our complex phase is just a incoherent mixture of the 
two ground states. (Ill) The comparison of the phase di- 
agrams of the renormalization group analysis and mean- 
field theory can be illustrated in Fig. (j4j. In order to 







QFIS 


Real BEC 2 




Complex BEC 




^^^^^ QFIS B 


NoBEC 






Real BEC 1 



-0.5 0.0 0.5 1.0 

Pi 

FIG. 4. (Color online) Comparison of phase boundaries from 
the renormalization group analysis and the mean-field anal- 
ysis. "MF" and "RG" indicate the boundaries from the 
mean-field and renormalization group analysis respectively. 
The complex BEC phase in renormalization group analysis is 
larger than the one from mean-field theory by a region named 
"QFISB". The coupling constants are g\ — §2 = 0.5, <?3 = 0.3 
and §4 — 0. 

compare the two phase diagrams with the same circum- 
stance we set <?4 = in the mean-field phase diagram. 
The expressions of the mean-field phase boundaries are 
b- = 2|a and §a = 4§- Fl. In Fig.ffl it's obvious to 

Ml gs Ml 2gi i-^J o 1—* 

see that the crucial difference is that the complex BEC 
phase is enlarged and the real BEC phase is suppressed 
in the RG phase diagram. 

In essence, the above differences can be explained as 
effects of the "quantum fluctuation induced symmetry 
breaking". The phase transition from the real BEC to 
complex BEC indicates the U(l) symmetry breakdown 
of field (j>i or <f>2- Mean-field description of the symme- 
try breaking is based on semiclassical approximation. In 
other words, it's a tree-level result. When we take into 
account the quantum fluctuation, the one-loop correc- 
tions can significantly change the model parameters and 
make some region of the real BEC phase become unsta- 
ble. In Fig. (@} these regions are labeled by "QFISB" . 
The possibility that the quantum fluctuations may pro- 
duce spontaneous symmetry breaking was first discussed 
by S. Coleman and E. Weinberg Q- They investigated 
a theory of a massless charged meson coupled to the elec- 
trodynamic field using effective potential method. Start- 
ing for a model without symmetry breaking at tree level 



they found that the one-loop effective potential indicated 
that a new energy minimum was developed away from the 
origin. This symmetry breakdown of the massless scalar 
field is completely induced by the quantum fluctuations 
of the electromagnetic field. Our theory reaches a qual- 
itatively analogous result using renormalization group 
analysis. One boson field plays the role of the electro- 
magnetic field and causes the symmetry breakdown of 
the other boson field through the interactions between 
them. 

Experimental proposal. — We finally consider the exper- 
imental feasibility. One leading candidate to observe 
QFISB in p-orbital BEC is 87 Rb atoms in a bipartite 
optical square lattice (ij. The optical potential can be 
constructed by crossing two laser beams with wavelength 
A = f,064nm and f/e 2 radius wq — fOO^m. The opti- 



cal potential reads — 



"0 ^[(.zcosa + y sma)e' lkx + 
eze —\ + e l ° z(e' lky + eer tk y)\ 2 . e < I and rj < I describe 
the imperfect reflection and transmission efficiencies, re- 
spectively. The typical values of e and r\ are e « 0.8f 
and r\ « 0.95. A BEC of 2 x I0 5 87 Rb atoms (in the 
F = 2, rn,F = 2 state) is produced in the optical trap. 
With Vo set to Vb/£rec = 6.2 the excitation of p-orbital 
band can be obtained by ramping 9 from 0.387T to 0.537T 
within 0.2ms. The initial values of chemical potentials 
Hi and \xi can be properly chosen by tuning a, which is 
the angle between the z axis and the linear polarization 
of the incident beam. 

We choose a = to prepare the system in the "QFISB" 
region of Fig. (QJ. The interaction strengths gi can 
be varied by changing the optical potential Vo [3, Hi- 
With different interaction strengths the momentum spec- 
tra can be obtained by taking the absorption images. 
For small interactions sharp peaks arise at momenta 
Kr +\ and K( + _^ in the Brillouin zone, where Kr±^±\ = 
1/2(±27T^/A,±27rfi/A) [|. This indicates a real BEC 
phase. For strong interactions sharp peaks should appear 
at the four momenta Kr+ t+ \, .KV Kr +i _\ and + ), 
which shows the appearance of the complex BEC phase. 
This would be a clear evidence of a phase transition from 
the real BEC to the complex BEC when we enhance the 
interactions. 

Conclusions. — In summary, we have investigated the 
quantum phase transition of the p-orbital boson gas in 
a two dimensional optical lattice using renormalization 
group method. We find that in contrast to the well- 
known "quantum triviality" shown by the s-orbital boson 
system the p-orbital one exhibits non-trivial behaviors. A 
phase transition from the real BEC phase to the complex 
BEC phase can be induced by the interactions between 
p x and p y orbital bosons. This non-trivial phenomenon 
is an effect of the quantum fluctuation induced symmetry 
breaking. Hence, the p-orbital boson condensates can be 
a counterexample of the quantum triviality exhibited by 
general two dimensional bosonic systems. 
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